A GENERALIZED BALANCED PAIR ALGORITHM 



BRIAN F. MARTENSEN 

Abstract. We present here a more general version of the balanced pair algorithm. 
This version works in the reducible case and terminates more often than the standard 
algorithm. We present examples to illustrate this point. Lastly, we discuss the features 
which lead to balanced pair algorithms not terminating and state several conjectures. 



1. Introduction 



The balanced pair algorithm was first introduced by Livshits |Liv-87| ILiv-92j for check- 



ing pure discrete spectrum for the Z-action of a substitution. A version of this algorithm 
was presented by Sirvent and Solomyak Sir-Sol-02j for irreducible substitutions. For sub- 
stitutions of Pisot type, Sirvent and Solomyak also give an explicit relationship between 
this algorithm and an overlap algorithm used in |Sol-97j for checking the pure discrete 



spectrum of the IR-action on the substitution tiling space. Recent results of Clark and 
Sadun ICla-Sad-TlS] give conditions for the conjugacy of the IR-actions on substitution 
tiling spaces which otherwise do not differ combinatorially or topologically. Their results 
give an immediate relation between the Z-action on the sequence space and the IR-action 
on the tiling space. The conditions given in Cla-Sad-'03] also allow us a way to generalize 



the balanced pair algorithm to reducible substitutions. The procedure for doing this is 
described in Section El below. 

The purpose of extending this algorithm two- fold. First, it creates a way in which 
the algorithm will terminate more often. When the balanced pair algorithm terminates 
has been much studied. Hollander |Hol-96j (see also )Hol-Sol-^T] ) has shown that the 
balanced pair algorithm will terminate for all 2-symbols Pisot substitutions. This fact 
along with the solution to the coincidence conjecture for 2-symbols [.Bar-Dia-02j . has 
shown that all 2-symbol Pisot substitutions have pure discrete spectrum. This has not 
yet been shown for an arbitrary number of symbols. In fact, it is not yet known whether 
the balanced pair algorithm terminates for all Pisot substitutions. 

The second reason for extending the algorithm to the reducible case is related to 
collaring or rewriting substitutions to obtain new (yet conjugate) substitutions. Collaring 
or rewriting procedures generally increase the number of symbols which turns irreducible 
substitutions into reducible ones. It would be beneficial to know that such procedures did 
not change the potential for a balanced pair algorithm to terminate. In particular, there 
are rewriting procedures which automatically produce coincidences. Thus the question 
of pure point spectrum (and even the coincidence conjecture) relies entirely on whether 
these reducible systems terminate. We explore these and other questions in Section |3] 
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2. Preliminaries 

Let A = {l,2,...,n} be a finite alphabet; A* will denote the collection of finite 
nonempty words with letters in A. A substitution is a map cp : A ^ A* and extends 
naturally to Lp : A* —>■ A* and : A^ A^ by concatenation. 

We associate to each word w a population vector 

which assigns to each Pi{w) the number of appearances of the letter i in the word w. 
To a substitution (p there is an associated transition matrix = {aij)i£A,j£A which 
ttij = Piifij))- Note that Ay,{p{w)) = p{(p{w)). We say a substitution tp is primitive if 
(p"^{i) contains j for all i,j & A and sufficiently large m. Equivalently, (p is primitive if 
and only if the matrix is eventually positive (there exists m such that the entries of 
A^ are strictly positive), in which case A^ has an eigenvalue A<^ larger in modulus than 
its remaining eigenvalues called the Perron- Frobenius eigenvalue of A^ (and ip). 

We form a space fl^ as the set of allowed bi-infinite words for tp. That is, u G fi<^ if and 
only if for each finite subword w of u, there are i & A and n G IN such that w is a subword 
of 'p"'{i)- We give the subspace topology (of A^ with the product topology) and denote 
the natural shift homeomorphism by a. Then {Q^f,,a) is a topological dynamical system 
which is minimal and uniquely ergodic provided the substitution is primitive. We will 
often use the fact that, due to unique ergodicity, any allowable word w appears in any 
u E fl^ with a well-defined and bounded positive frequency. 

The spectral type of the measure-preserving transformation [Q^p, cr, //) is, by definition, 
the spectral type of the unitary operator U^, : /(■) t-^ /(o"-) on L'^{fl^,fi). We say 
(n^, a, fj,) has pure discrete spectrum if and only if there is a basis for L'^^Q^, /i) consisting 
of eigenfunctions for U^,. 

Additionally, to each substitution ip and any u G fi<^ we can form a tiling of the real 
line by intervals. For each letter i E A we assign a closed interval of length We refer 
to the vector L = (/i, . . . , /„) as the length vector. We then "lay" copies of these closed 
intervals down on the real line (so that they do not overlap on their interiors) according 
to the order prescribed by u = . . . U-i.uqUi . . ., with the placement of the origin given by 
the decimal point. There is a natural translation acting P* on a tiling T which forms new 
tilings by simply moving the origin by a distance t G IR to the right. For a fixed u, we 
form a compact metric space by taking the completion of the space of the translates of u. 
Here the completion is taken with respect to the metric which considers two tilings to be 
close if they agree on a large neighborhood about the origin after a small translation. If 
the substitution is primitive, then this space is independent of u, is minimal and uniquely 
ergodic. We refer to the H-action on (n;^,P*,L) as the tiling dynamical system. It is 
topologically conjugate to the suspension flow over the Z-action (fi^, a), with the height 
function equal to k over the cylinder i. The system (fi^, P*, L, /i) is then said to have 
pure point spectrum if there is a basis of L^(f2(^, P*, L, fi) consisting of eigenfunctions for 
the IR-action. 

A primitive substitution is of Pisot type if all of its non-Perron-Frobenius eigenvalues 
are strictly between and 1 in magnitude. Solomyak and Sirvent |Sir-Sol-02] have shown. 



for a substitution of Pisot type, that if the IR-action has pure discrete spectrum then so 
does the Z-action. A recent result of Clark and Sadun |Cla-Sad-03j implies that this can 
be strengthened to an "if and only if statement. 
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In fact, the result of Clark and Sadun has more general applications than just the Pisot 
case and we briefly describe some aspects of their results. First, however, we should give 
a few observations in order that we may put their results in context. 

We denote as L;^ the left eigenvector associated to the Perron-Frobenius eigenvalue A 
of the transition matrix A^. The choice of as a length vector causes any tiling which 
is (combinatorially) fixed under ip to be self-similar. The similarity map is expansion 
by A and the tile-substitution is the geometric version of (p. On the other hand, the 
length vector Li = (1, 1, . . . , 1) is a natural choice since the group action of the time-1 
return map for the tiling dynamical system is conjugate to the shift homeomorphism 
of the substitution dynamical system. Other length changes, though producing tiling 
dynamical systems conjugate to the suspension of the shift, will produce no such group 
action. It is easy to see that the IR-action of the tiling dynamical system with constant 
length vector Li has pure discrete spectrum if and only if the Z-action does. Therefore, 
questions about the relation of the H-action with length vector L to the Z-action can 
be rephrased into questions about the conjugacy of the IR-actions of the systems using 
L and Li, respectively. 

Clark and Sadun |Cla-Sad-'03] give explicit conditions for the IR-action associated to 
one length vector to be conjugate (up to an overall rescaling) to that of another via a 
homeomorphism which preserves the combinatorics (i.e. is homotopic to the identity). 
Assume the length vectors L and L' have been rescaled to agree in the Perron-Frobenius 
direction. Then, the conjugacy occurs if and only if hA^ — h'A^ ^ as n —* oo. Note 
that every length vector can be written as a linear combination of vectors living in the 
Perron-Frobenius eigenspace, the small eigenspaces (those associated to the eigenvalues 
of magnitude less than one) and the large eigenspaces (those associated to eigenvalues 
of magnitude greater than or equal to one, though we don't include the P. F. -eigenspace 
here) . In particular, after the rescaling, L — L' must avoid the large eigenspaces for the 
conjugacy to occur. If the substitution is of Pisot-type, then all systems length changes 
produce conjugate systems and we arrive at the result above. 

We point out a few special kinds of substitutions we will be considering. A substitution 
ip is said to be irreducible if the characteristic polynomial of its transition matrix is 
irreducible. In particular. La of an irreducible substitutions is such that none of its 
entries are rationally related to one another. Lastly, a substitution ip is said to have 
constant length if the number of letters appearing in (p{i) is the same for all i E A. 



3. Balanced Pair Algorithms 

We now describe the balanced pair algorithm in a variety of circumstances. We begin 
with the case that the substitution is irreducible. This case was studied extensively in 
|Hir-.^ol-n2| as an adaptation of the algorithm of |Liv-92j . The extension of this algo- 
rithm to a specific class of reducible substitutions containing "letter equivalences" can 
then easily be described. We end this section with description of how one extends this 
procedure for generic reducible cases. 



3.1. Irreducible case: A pair of allowable words u and v are called balanced if the have 
the same population vector. We write I " I if and v are balanced. Note that if I " I is 

balanced, then so is "^["^ . 
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Let the right infinite sequence u = ^0^1 • • • be fixed under the substitution and let 
w he Si non-empty prefix of u. Since w appears in u with positive frequency, u can be 



written as u = 10X110X2 
in the following way: 



We may then speak of splitting 



into balanced pairs 



u 




wXi 




WX2 






Xiw 




X2W 



Since appearances of i« in u are bounded, there are only finitely many different balanced 
pairs encountered in the process above. We may further reduce each of these to form a 
finite set of irreducible balanced pairs which we will refer to as Ii{w). 
We may now inductively define, for n > 1: 



w 



{ 


u 




u 




V 




V 



in the reduction of 



appears as an irreducible balanced pair 



for some 



y 



el 



Let I{w) = U^i Iniw)- If liw) is finite, then we say that the balanced pair algorithm 
associated to a prefix or bpa-t/;, terminates. Below, we will state how this algorithm is 
used to determine pure discrete spectrum, though our main interest here is in determining 
precisely when the algorithm terminates. We illustrate the computation of I{w) for a 
simple example. 

Example 3.1. 

Consider the substitution given by: 



Then u 



Now, 



1 
2 

11211212112 ... and we take w = l. 



112 
12 

Then it is easy to see that: 
12 
21 



1 




1 




2 


and 


12 




1 




12 




1 




2 






























1 




1 




2 


21 




1 




21 




1 




2 



so that I2 ■ 
terminates. 



{i;i 



1 12 1 12 
21 1 ' 1 2 



} and further, I2 = I3 



I{w). Thus, the algorithm 



3.2. Substitutions involving letter equivalences: The balanced pair algorithm as 
originally described by Livshits |Liv-87t ILiv-92j has an additional feature. Consider two 
letters i,j E A to be equivalent if and (p'^{j) have the same number or symbols 

for all n G IN. Then, a pair of words is balanced if both contain the same number of 
symbols from each equivalence class. The algorithm runs in the usual way, starting with 
an initial list of balanced pairs, substituting and reducing. It stops if no new balanced 
pairs are produced. In particular, in the constant length case all letters are equivalent 
so all irreducible balanced pairs are just pairs of symbols. Thus, the algorithm always 
terminates in this case. This way, the algorithm includes Dekking's criterion pek-78j in 
the constant length case. For any irreducible case, no two symbols are equivalent and 
thus the algorithm runs just as before. 
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Example 3.2. (Substitution of constant length) 
Consider the substitution given by: 

1 ^ 112 

2 ^ 122 

which is a constant length substitution. If we ignore the equivalence relation of Livshits, 
then the troubling balanced pair is \ ^\. Iterating this pair, we see: 



12 




1 




1212 




2 


21 




1 




2211 




2 



1212 




1 




1212211212 




2 


2211 




1 




2212211211 




2 



In particular, this process generates new balanced pairs of the form | | for longer and 
longer words w. 

Once we take the equivalence relation into account, however, all letters are equivalent 
and therefore the above process terminates. 

Example 3.3. (A nan- constant length substitution) 

Consider the substitution given by: 

1 ^ 31 

2 ^ 412 

3 ^ 312 

4 ^ 412 

If we were to again ignore the equivalence relation of Livshits, then there is a potential 

problem with balanced pairs that "match up" the letter 3 with 4. To see this, note that 
the right eigenvectors of the transition matrix in some sense describe the frequencies in 
which letters appear. Here, there is one large eigenvalue other than the Perron-Frobenius 
eigenvalue. The eigenvalue is 1 and it has right eigenvector (0,0, 1,-1)^. We can thus 
view a balanced pair of the form | ^" [ as initially having an abundance of a 3 and a lack 
of a 4 on top. Since this corresponds to the right eigenvector of 1, this difference persists 
under substitution. Note that we say this is a "potential" problem as this association 
to a large eigenvalue in itself will not force the bpa-w to not terminate. In this case, 
however, we show that this does in fact occur. After shifting the fixed word u = 312 . . . 
with w = 31, we see the balanced pair I ^\t^ |. Iterating this pair, we see: 



31412 




41231 





3123141231412 
4123141231231 



Thus, we generate new balanced pairs of the form | 4^231 1 ^'^^ longer and longer words w. 
The equivalence relation tells us however that 2, 3 and 4 are actually equivalent letters 
and thus the above process terminates as | ^ | and | 2 1 1 3^ | are balanced. Note that the 
left Perron-Probenius eigenvector for this system is (1, A, A, A) so that 2, 3 and 4 all have 
the same length in the self-similar system. 
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3.3. Reducible case: The example above gives some hint as to how to describe the 
balanced pair algorithm in the generic reducible case. For a system whose length of tiles 
have rational relations which persist under substitution, it may be the case that entire 
words should be identified even when no individual letters are. We therefore introduce 
an equivalence relation which exploits these rational relations. Consider two words v and 
w to be equivalent if: 

L(p(<^"(^;)) - p(V9"H)) = 0, Vn e IN, 

or equivalently, if 

L{A;piv)-A;piw)) = 0,\/nef^. 

We write v ~l w to emphasize the dependence of the equivalence relation on choice of L. 
In the case in which we use the left Perron- Frobenius eigenvector for our lengths, the 
equivalent pairs correspond simply to the geometric balanced pairs of |Sir-Sol-02] . In fact, 
we will show below that in this case, the algorithm is equivalent to the overlap algorithm 
and for a general L is meant to bridge the gap between the balanced pair and overlap 
algorithms in the reducible case. 

The algorithm once again runs in the usual way, with a pair of words considered 
balanced if they are equivalent under our relationship above. We denote this algorithm 
by bpa{w,L) and the set of balanced (or equivalent) pairs by I{w,Li) to emphasize that 
their is now an additional dependence on L. 

Remark 3.4. Our equivalence condition above that Ij{p{ip'^{v)) — p{f"'{w))) = 0, Vn G IM, 
could be weakened to allow L{p{Lp^ (v)) — p{Lp"-{w))) — as n ^ oo, which essentially 
allows inconsequential length changes in the small eigenspaces mentioned above. 

Example 3.5. 

Consider the substitution given by: 

1 ^ 112 

2 2321 

3 12. 

The eigenvalues of the transition matrix are ^^^^ and 1 while = (1, ~^+^^ ^ ^~^^ ). 
The right eigenvector associated to the eigenvalue 1 is (2, —1, —1)"^. This indicates that 
there is a potential problem with balanced pairs which "match up" the word 11 with 
23. Such a matching occurs in the balanced pair | ^32^^ |, which was found by performing 
the balanced pair algorithm after shifting the fixed word by two spaces. Iterating this 
balanced pair, we generate balanced pairs of the form | 23^^'^ | for longer and longer words 
w. 

Fortunately, La ■ (2, —1, —1) = so that 11 and 23 are equivalent words and hence | 23 | 
is balanced and the above algorithm will terminate. Also note that Li is also perpendic- 
ular to (2,-1,-1) so that the algorithm will also terminate for what we consider to be 
the other interesting case. This is immediate from the fact that these two systems are 
conjugate by |Cla-Sad-'03] . On the other hand, L = (1, 1, 2) for example is not conjugate 
to these and furthermore the algorithm does not terminate in this case. 

Proposition 3.6. Let S{L) = {{u,v) : u ~l v}. Then, for any length vector L, S{L) C 
^(La) 
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Proof. Let u ~l v. Let z = p{u) — p{v). Since L has only positive terms, it must have 
a component in the La direction. Let L = axLx + Yl]=i ^j^j + J2j=i '^i'^ji "where Bj is 
a left eigenvector for eigenvalue j3j > 1 and Sj is a left eigenvector for eigenvalue aj < 1 
for each j. Then, 

I s 

If La • z 7^ 0, then the above implies A" ~ X]j=i ^jPj' large enough n and some 
constants Cj. This is a contradiction since the Perron- Frobenius eigenvalue dominates 
all others. □ 

Corollary 3.7. For any length vector L, if the bpa{w,L) terminates, then the bpa^WjLx) 
terminates. 

Let C denote the span of the right eigenvectors with eigenvalues greater than or equal 
to 1 in magnitude but strictly less than the Perron- Frobenius eigenvalue. Similarly, 
let S denote the span of the right eigenvectors with eigenvalues strictly less than 1 in 
magnitude. We say that {u,v) lies in a vector space V if {p{u) — p{v)) G V. 

Remark 3.8. Then, the set of equivalence words lies entirely in the spaces since they 
form precisely what is perpendicular to La. In our examples above, choosing an L which 
missed identifying equivalent pairs C would have led to the bpa(i«,L) not terminating. 
Contrast this with Example 13.31 in which choosing an L which neglects to identify 2 and 
4 will have no effect on whether the algorithm terminates. The vector associated to the 
pairing of 2 and 4, namely (0, 1, 0, —1) lives in the zero-eigenspace and thus differences 
in frequencies should be quickly dispelled. Our contention here is that generally the 
equivalence relations which live in S should not affect the algorithm whereas those in £ 
affect it greatly. 

Corollary 3.9. u ~l v implies p{u) — p{v) G £ © 5. The converse is true z/L = La. 

A balanced (equivalent) pair |*|, for i G is called a coincidence. We say that a 
balanced pair | " | leads to a coincidence if there exists m such that the reduction of 

contains a coincidence. Notice that coincidences lead to coincidences since 



las nothing but coincidences in its reduction 



Theorem 3.10. Let ip be a primitive substitution such that u = UqUi . . . is a right infinite 
fixed word and let L be a length vector. 

(a) If for some prefix w the bpa{uj, L) terminates and every equivalent pair in I{vu, L) 
leads to a coincidence, then {Q^,r^ ,Ijx) has pure discrete spectrum. 

(b) // the bpa{w, L) terminates for some prefix w = uq . . . Um such that Um+i = uo, 
and (f2^,r*,LA) has pure discrete spectrum, then every balanced pair in I{w,L) 
leads to a coincidence. 

Before beginning the proof of Theorem I3.1UI we make the following observations re- 
garding the densities of coincident pairs. Let u = uqUi ... be a fixed right-sided sequence. 
Let z he a prefix of u and u — L ■ p(2;) = vqVi . . .. Let D{z) = {ui : Ui = vj some j 
with uo...Ui^i ~L vo...Vj-i} Suppose we define a density function, densi,{D{z)) = 
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limj.^;^ ^•p(-PW'~'"o---ffc) ^Yie limit exists. The existence of this hmit follows from the 

unique ergodicity of F*, L). Notice that for L = Li, this definition of density agrees 
with that used in |Sir-Sol-02] for the (irreducible) balanced pair algorithm and for L = La 



it agrees with that of |Sol-97j for the overlap algorithm. Now, by Proposition 13.61 coin- 
cident pairs for S(L) are also coincident pairs for £{Lx). 

Remark 3.11. The proof of this theorem differs from that of the irreducible case only in 
the way in which we define density and the set of irreducible balanced pairs. We therefore 
only include a sketch of the proof that follows closely a sketch provided in |Sir-Sol-0'2] 



for the irreducible case. The full details of that case has been worked out in |Hol-96j . A 
theorem of this sort was proved by |Liv-92j . though coincidences and the balanced pair 
algorithm go back to jDek-78j and |Mic-78j , respectively. Part (a) is largely contained in 
|Que-87| . 

Proof, (of Theorem I3.1U|) 

Let w be a prefix of the fixed word u. Denote by D{z) the density defined above. We 
will be interested in densi,{D{ip''{w))) as / ^ oo. Let 

be the reduction of | ^p^^ | into irreducible equivalent pairs. For an equivalent pair /3 = | " | 
let |/3| = L ■ p(m) and (5(/3) = {ui E u : Ui = Vj some j with uo . . . ~l vq . . . vj^i}. 
Then 

dens^iDi^ (w))) = Jim _ ■ 

Now consider the substitution on the set of irreducible balanced pairs I{w,Ij). By 
definition clearly u*^') G I{w, L)^ and u*^') = (p''(n^^^), where u^^^ is the reduction of | ^|"|^| 
into irreducible equivalent pairs. 

There is a directed graph Q{ip) associated with the substitution 0. Its vertices are 
labelled by the members of I{w, L), and for every vertex (3 there are directed edges from 
f3 into of the letters of 0{l3) with multiplicities. 

(a) Let (3 G I{w, L) be an irreducible equivalent pair which is not a coincidence. By 
assumption there is a path in the graph Q{0) leading from /3 to a coincidence. Since all 
the edges from coincidences lead to coincidences, a standard argument shows that the 
frequency of the symbol /3 in u*^'^ = (p''{u^^^) goes to zero geometrically fast as / — oo. 
Since I{w, L) is finite and 1 — 6{l3) > if and only if /3 is a non-coincidence, it follows 
that 



1 — dens{Dp^) < const ■ 7' 

for some 7 G (0, 1). This implies that u is mean-almost periodic so we can conclude (sim- 
ilar to |Que-87| ), VI.25) that ip has pure discrete spectrum. We note that the argument 
only relied on I{w, L) being finite and that the equivalence relation does not create false 
coincidences. 

(b) Let w = Uo-.-Um be such that Um+i = Uq and pi = h ■ p{ip\w)). It follows 
from |Hos-8f)j (see also |Sol-97j . Theorem 4.3) that lim^^oo e^^^'^^' = 1 for any eigenvalue 
^2mx q£ ^Yie dynamical system (i^;^, F*, L, /i). If the spectrum is pure discrete, then the 
eigenf unctions span a dense subset of L^(fi<^) so that lim^^^oo \ \U^''f — /II2 = for every 



GENERALIZED BALANCED PAIR ALGORITHM 



9 



/ G Taking / to be the characteristic function of the cyhnder set corresponding 

to z G ^ with heights prescribed by L, we obtain (see Que-87| , just after Lemma VI. 26) 



limi^oo densi,{D{Lp\w))) = 1. On the other hand, suppose that there is an irreducible 
equivalent pair in I{w, L) which does not lead to a coincidence. Then there exists an 
irreducible component Qo of the graph Q (0) which contains no coincidences. There exists 
/o such that for every / > Iq elements of the component Qq occur in u^'^ with positive 
frequency. (Note that different elements of Qq may occur for different /). Further, it can 
be shown that this frequency is bounded away from zero as / — > oo. Since 1 — 6{/3) > 
for all /3 G Qo, it follows that dens-L,{D{ip'- (w))) -/-^ 1, which is a contradiction. 

□ 

We also note the following relationship between the balanced pair algorithm and the 
overlap algorithm in the case that the tiling space is self-similar. 

Proposition 3.12. Let if he a 'primitive substitution such that </?(!) begins with 1. Then 
the overlap- algorithm associated to x = x{vu) terminates with half- coincidences if and 
only if bpa{w, La) terminates. 

Proof. Assume the bpa(w, L^) terminates. Then the distance between half-coincidences 
(endpoints of our equivalent pairs) arising from looking at (T, T — A"x) is bounded, where 
X = x{w) = L • p{uj). Theorem 5.6 of |Sir-Sol-02] implies that the overlap algorithm 
associated to x{w) terminates with half-coincidences. 

Assume the overlap-algorithm associated to x = L-p(iy) terminates with half-coincidences. 
Again by Theorem 5.6 of [ Sir-Sol-02j . the distance between half-coincidences arising from 
(T, T — X^x) is bounded. Hence bpa(w, L^) terminates. 

□ 

Example 3.13. A non-terminating example 

We now give an example which will not terminate even with the extended equivalence 
relations presented here. This example is a rewriting of the Morse-Thue systems in which 
we have forced coincidences. We will use La as our length vector so that we are conjugate 
to the original Morse-Thue system and hence do not have pure discrete spectrum. This 
example therefore cannot terminate for any version of the balanced pair algorithm. 

The substitution is given by: 

1 1234 

2 ^ 124 

3 13234 

4 1324 

Here, the P.F. eigenvalue is 4 with left eigenvector (3, 2, 4, 3). (Note that this vector also 
corresponds to the Z-action of the original Morse-Thue system.) The other important 
eigenvalue (the remaining two are zero) is 1 with right eigenvector (1, 1, —2, 1)-^. Thus, 
abd is equivalent to cc, however these words do not cluster close enough to each other to 
aid in terminating the balanced pair algorithm. Notice also the system generated by the 
length vector (1, 1, 1, 1) is not conjugate to the system generated by La. But by Corollary 
13. 7| the balanced pair algorithm will also not terminate for this length vector. In fact, 
using techniques from |Cla-Sad-'03] to directly compute its spectrum, one can see that it 
will not have pure discrete spectrum either. 
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4. Open problems and conjectures 

The motivation for studying the effects of the new equivalence relation on the balanced 
pair algorithm was not only to understand reducible substitutions, but was mainly to aid 
in analyzing substitutions which have been collared or rewritten. We use the example 
below as a test case for a general Pisot substitution. Beginning with a Pisot substitution, 
we will rewrite it into an equivalent substitution that always begins with the same letter. 
In this way we force every balanced pair to lead to a coincidence. It remains to show that 
the new substitution will terminate in order to show that all Pisot substitutions have pure 
discrete spectrum. A difficulty arises in that rewriting increases the size of the alphabet 
and can therefore add additional eigenvalues of and ±1. The zeros do not concern us, 
but the roots of unity might. Considering the Z-action on our original substitution by 
changing tiles to unit lengths produces an integer vector in the rewritten substitution 
which will not have a component in the eigenspace of the roots of unity. Further, these 
eigenspaces generate equivalent words so that it is our hope that the equivalence classes 
will always force the algorithm to terminate. 

Example 4.1. A rewritten Pisot substitution 

Consider the substitution (f given by a ^ abb and b ba. Using the rewriting 
procedure of |Bar-Dia-'OT] . we square (p and generate a substitution on 1 = abbb,2 = 



ab, 3 = aabbb and 4 = aabb by: 

1 122334 

2 1224 

3 12322334 

4 1232234 

The eigenvalues of the transition matrix of this substitution are 0, 1 ± ^ and 1. omce 
the original substitution is Pisot, it is insensitive to length changes in a and b. The 
new substitution will also be insensitive to any length change which is consistent with 
length changes in a and b. For example, setting a = b = 1 generates the length vector 
L = (3, 2, 5, 4) and will in fact miss the eigenspace associated to 1. (Note, however, that 
Li = (1, 1, 1, 1) does not.) The right eigenvalue of 1, Vi = (1, 1, —2, 1)-^, is perpendicular 
to L so that abd ~l cc. Because Li is not perpendicular to Vi, this makes the hpa.{w, L) 
difficult to run. Some tedious calculations reveal that hpa.{w, L), where one must keep an 
eye out for the equivalence relations, will terminate. This example produces 30 different 
irreducible balanced pairs, the longest one of which contains a word of length 11. We 
suspect that bpa(ty,Li) will not terminate, but this has also not yet been shown. 

Generally, if the balanced pair algorithm terminates for a Pisot substitution, |Sir-Sol-02j 
gives us that the distance between half-coincidences is bounded. For the rewritten sub- 
stitution, the new tiles are compositions of smaller old tiles and these half-coincidences 
may occur "internally." We suspect this cannot happen and that the algorithm must 
terminate for the new substitution as well. More precisely we have: 

Conjecture 4.2. Let (f be a Pisot substitution on n-letters. Let (p be a rewriting of ip 
so that or some letters b, e in the rewritten alphabet, (p{i) = b . . .e for all i. Assume that 
the balanced pair algorithm for the original Pisot substitution ip terminates. Then the 
balanced pair algorithm for the rewritten substitution (p also terminates for length vector 
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An immediate corollary of this would be that if the balanced pair algorithm of a Pisot 
type substitution terminates, then it must do so with coincidences. 
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